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ABSTRACT 

Collimated fast winds (CFWs) have been proposed to operate during the post-AGB 
evolutionary phase (and even earher during the late AGB phase), responsible for the 
shaping of pre-planetary nebulae (PPNs) and young planetary nebulae (PNs). This 
paper is a follow-up to our previous study of CFW models for the well-studied PPN 
CRL 618. Previously, we compared our CFW models with optical observations of CRL 
618 in atomic and ionic lines and found that a CFW with a small opening angle can 
readily reproduce the highly collimated shape of the northwestern (Wl) lobe of CRL 618 
and the bow-like structure seen at its tip. In this paper, we compare our CFW models 
with recent observations of CRL 618 in CO J=2-l, J=6-5, and H2 1-0 S(l). In our 
models, limb-brightened shell structures are seen in CO and H2 at low velocity arising 
from the shocked AGB wind in the shell, and can be identified as the low-velocity 
(LV) components in the observations. However, the shell structure in CO J=2-l is 
significantly less extended than that seen in the observations. None of our models can 
properly reproduce the observed high- velocity (HV) molecular emission near the source 
along the body of the lobe. In order to reproduce the HV molecular emission in CRL 
618, the CFW is required to have a different structure. One possible CFW structure 
is the cylindrical jet, with the fast wind material confined to a small cross section and 
collimated to the same direction along the outflow axis. 

Subject headings: planetary nebulae: general - stars: AGB and post-AGB — stars: 
mass loss - stars: winds and outflows 



1. Introduction 

At the end of the evolution of low- and intermediate-mass stars, pre-planetary nebulae (PPNs) 
are seen associated with post asymptotic giant branch (post-AGB) stars. They will turn into 
planetary nebulae (PNs) in less than 1000 years as the post-AGB stars become hot white dwarfs. 
Their shaping mechanism is still unclear and is closely related to the mass-loss processes during 
the end phases of the evolution. 
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Many PPNs and young PNs were found to have highly aspherical shapes, with a significant frac- 



tion h aving highly coUimated bipolar or mu l tipolar lobes (ICorradi fc: Schwardll995l : ISchwarz et. al 



19971 : ISahai fc Trauger I Il998l : ISahai 



2001 



Sahai et al 



200 



Point symmetry, rather than ax- 
isymmetry, b etter characterizes the g eometry of the majority of these objects [striking exampl es 
are He 2-115 fISahai fc Trauger 1 1 19981 ) and He 3-1475 teorkowski. Blondin. fc Harrington ! 1 199?! )]. 
As a result, instead of sp herical fast winds as in the generalized interacting stellar winds (GISW) 
model (see e.g., review by lBalick fc Frank II2002I ). colhmated fast winds have been proposed to op- 
erate during the post-AGB phase (and even earli er during; the late AGB phase) , and b e the primary 



agents for the shaping of PPNs and young PNs (jSahai fc Trauger Ill998l : ISahai Il200ll ). Collimated 



fast wind (CFW) models have been used to account for the morphology and kinematics of a few 
well-studied PPNs a nd PNs, with some assuniing: a radial wind with a small opening angle (Lee & 
Sahai 20 03. hereafter [ Paper AlAkashi fc Sokerll2008l). some assuming a cylindrical jet either unmag- 
netized (Iciiffe et al.lll995l : ISteffen fc L6pezlll998l : iGuerrero et al.ll2008l ) or mag netized (ILee fc Sahai 



2004), and some assuming a bullet (a massive clump) along the outflow axis ( Dennis et al. 20081 ) 



This pap er is a follow-up to our previous study of CF W models for th e well-studied PPN CRL 



618 (jPaper j ). CRL 618 is located at a dista nce of 900 pc (Goodrich 



lobes at different orientations in HST imag es (ITrammel" 



tipolar" morphological classification (Sahai et al.l 120071 ). perhaps resulting from multiple ejections 



fc Goodrich 



1991). It shows several narrow 
2002 ) and belongs to the "mul- 



at different orientations. The general structures of the different lobes are similar and we f ocus only 
on the northwestern (Wl) lobe, which seems to be better separated from other lobes. In IPaper I . 
we compared our CFW models with optical observations of CRL 618 in atomic and ionic lines and 
found that a CFW with a small opening angle can readil y reproduce the hig:hly collimated shape 
of th e Wl lobe and the bow-like structure seen at its tip (ISanchez Contreras. Sahai. fc Gil De Paz 
2002). However, it may have difficulties in reproducing properly the high- velocity optical emission 
along the body of the lobe. In this paper, w e compare our CFW models with recent ob servations 



of CRL 618 in H^ 1-0 Sfl ) (ICox et al 



2003), CO J=2-l dSanchez Contreras et al.ll2004l ) and J=6 



5 (INakashima et al.l 120071 ) . We find that our CFW models also have difficulties in reproducing 
the high-velocity molecular emission in the Wl lobe and the CFW is required to have a different 
physical structure. 



2. Numerical Methods and Physical Settings 



The two -dimensional hydrodynamic code, ZEUS 2D, is used for the simulations of our models, 
as in lPaper"^ . This code has been now modified to include molecula r cooling; and the tin ie dependent 
chemistry of hydrogen by solving the following equations (see also lSuttner et al.lll997l ): 
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where e, v, and T are the internal energy density, thermal pressure, velocity, and tempera- 
ture, respectively. Also, n is the hydrogen nuclei density, / is the fractional part of hydrogen 
molecules with riHs = (i-e., / = for atomic/ionic gas and / = 0.5 for molecular gas), and g 
is the fractional part of ionized hydrogen with nfj+ = gn. Helium is included as a neutral com- 
ponent with riHe — O.ln and thus n — p/(1.4mH), where p and mn are the mass density and 
the mass of atomic hydrogen, respective ly. Here D and R a re dissociation and ressociation rates 
of molecular hydrogen, respectively (see ISuttner et al.l 119971 . and reference therein), / and C are 
ionization and recombination rates of atom ic hydro gen, respectively. is the optically thin ra- 



diative cooling!; from a toms and ions (as in iPaper I ), wi th the cooling; at hig h temp erature from 
MacDonald fc Baileyl (Il98ll ) and low temperature from IPalgarno fc McCravl fll972l). A a/t is the 



optically thin radia t ive cooling; from molecules, including; H2 flHollenbach fc McKeelll979l ) and CO 
f|McKee et al.lll982l : iHollenbach fc McKeel [l989l ) . The chemistry of CO is not included. The CO 
abun dance is assumed to be constan t and equal to 2xl0~^ of the number of hydrogen molecules. 



as m 



Sanchez Contreras et al.l (120041 ) . We believe this assumption does not significantly affect our 



conclusions, it only affects the small-scale structure o f the CO emission. The equations of state 
and material coefficients are (for similar derivation see lSuttner et al.l 119971 ): 
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where fc, p, and Cy are the Boltzmann constant, mean molecular weight, and specific heat, respec- 
tively. A scalar color tracer c is also included in the simulations to track the fast wind, it is one 
for the fast wind, zero for the A GB wind , and a value between one and zero for a mixture of the 
fast wind and AGB wind. As in IPaper U the simulations are performed in spherical coordinates 
but presented in cylindrical coordinates (z, R), with the z-axis being the outflow axis. 



Our CFW models are based on models 1 and 4 in IPaper U which were found to be the best 
models for the Wl lobe of CRL 618 as seen at optical wavelengths. In these models, a CFW 
with a small opening angle is assumed to emanate radi ally froni the vicinity of the post- AGB star, 
interacting with a spherical AGB wind (for details see IPaper u ) . In our simulations, the CFW is 
launched from the inner boundary of the simulation domain, which is at 5x10^^ cm (^ 333 AU) 
away from the post-AGB star. The AGB wind is assumed to be molecular with a temperature of 
10 K and have a mass-loss rate of 3 x 10~^ Mq yr~^ with an expansion velocity of 20 km s~^. The 
CFW is assumed to have a mass-loss rate of 2.5x10"^ Mq yr~^, a speed of 300 km s~^, and an 
opening angle of 10°. It can either be steady or pulsed with a temporal variation in the density 
and velocity (see Table [1]). It can be either atomic or molecular, depending on the temperature of 
the CFW to be assumed. 
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3. Recent Observations in CO and H2 



CRL 618 has been recently observed in H2 1-0 S(l) at ^ 0(^5 resolution f|Cox et alJl2003l ) and 
in CO J=2-l (jSanctiez Contreras et al. 2004 ) and J=6-5 ( Nakashima et al. 200?! ) at ^ resolution 
(see Figs. [T]and[2j). It is multipolar, but here we focus only on its Wl lobe that extends ^ 6^' to 
the west from the source. The systemic velocity in this region is —21.5 km s~^ LSR and the Wl 
lobe is mainly redshifted with a velocity ranging from -40 to 150 km s-^ LSR. The CO and H2 
emission toward the Wl lobe can be roughly separated into two components: a slow or low-velocity 
(LV) component with a velocity lower than ~ 20—30 km s~^ from the systemic and a fast or high- 
velocity (HV) component with a velocity extending to ^ 170 km s~^ from the systemic. The LV 
component is extended. In CO J=2-l, it forms a limb-brightened shell structure around the optical 
lobe, extending to ^ 5^' away from the source. In CO J=6-5, it also forms a shell structure but 
on ly in the southern part of the lobe, extending to only ^ 3^' away from the source. As discussed 



m 



Nakashima et al 



20071 ). however, it should be more extended, because most of the extended 
flux has been resolved out in their observations. In H2, it is spatially unresolved, extending mainly 
from 2^' to &^ from the source, slightly ahead of that seen in CO J=2-l. On the other hand, the 
HV component is more compact. In CO J=2-l, it extends to ^ 2(^5 away from the source with 
the velocity increasing linearly with the distance from near the systemic velocity to the highest 
velocity. In H2, it is seen with three localized emission peaks separated by ~ 2'^ (1) a peak close 
to the source with a range of (blueshifted and redshifted) velocities, (2) a bright peak at ^ 2(^5, 
where the tip of the CO HV component is, with a broad range of velocities extending from near 
the systemic velocity to the highest velocity, and (3) a faint peak at ^ 4^^ at ^ 80 km s~^ LSR. 



4. Results 

In the following, we present our models and the comparison with the observations. Figure 
[3] shows the distributions of hydrogen nuclei density and temperature with molecular fraction in 
our models at the age of ^ 160 yrs, when the outflow lobe has a length of ^ 6000 AU, similar to 
the deprojected length of the Wl lobe of CRL 618. The separation between the AGB wind and 
the CFW is delineated by the color tracer c = 0.5, at which half is the AGB wind and half is the 
CFW. In order to compare with the observations, we also derive intensity maps of the LV and 
HV components (Fig. [4j) and position velocity (PV) diagrams (Fig. [5j) for the CO J=2-l, J=6-5, 
and H2 1-0 S(l) emission from our models, assuming the latter to be optically thin and arising 
from gas in local thermal equilibrium ( LTE). We assume a distance of 900 pc and an inclination 



of 30°, values appropriate for CRL 618 (ISanchez Contreras et al.ll2004l ). The HV components are 



derived by integrating the emission with velocity higher than 50 km s~^ from the systemic and 
the LV components by integrating the emission with velocity within 50 km s~^ from the systemic. 
Two angular resolutions, 0(^1 and 0(^5, are assumed, with the latter for comparing with the current 
observations. 
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4.1. Model 1: Steady Atomic CFW at 10,000 K 



Model 1 here corresponds to model 1 in lPaper II . The CFW has a temperature of 10 K and is 



thus atomic. As it blows into the AGB wind, it produces a collimated outflow lobe, which is a thin 
shell with a cavity (Fig. [3^). As shown, the shell consists of shocked AGB wind in the outer shell 
and shocked fast wind in the inner shell (Fig. [3]A.). Since the shell is already radiative and thus 
momentum-driven even without molecular cooling, additional molecular co oling; doe s not change 



the shell dynamics, it only reduces the shell thickness as compared to that in lPaper J . In the shell, 
the shock becomes stronger going from near the source toward the tip of the lobe. As a result, the 
temperature of the shock AGB wind increases from tens K near the source to above 10^ K near 
the tip (Fig. [Hb). In addition, the shocked AGB wind is mainly molecular near the source but 
becomes mainly atomic near the tip (except for the newly shocked AGB wind) because of shock 
dissociation. In contrast, the shocked fast wind is mainly atomic. 

In the simulation, a jet-like structure is seen at th e tip bec ause of the accumulation of material 
there due to the shock focusing effect, as discussed in IPaper^ . Such accumulation might be artifi- 
cially enhanced in our 2D simulations because material cannot flow across the symmetry axis. This 
enhancement in turn would lead to more cooling and further accumulation, and then an increase 
in the formation rate of molecules. Therefore, the molecular emission from the lobe tip might be 
highly overestimated and should be kept in mind when comparing our models to the observations. 

In this model, the LV components of the CO J=2-l, J=6-5, and H2 emission are seen forming 
limb-brightened shell structures around the cavity (Fig. [H model 1), as in the CO observations 
(Fig. [1]). They arise from the shocked AGB wind in the shell. Since the temperature in the shell 
increases with the distance from the source, the shell structure extends further and further away 
from CO J=2-l to CO J=6-5 and to H2. The CO J=2-l emission traces the cold (20-100 K) 
gas extending to ^ "2!' from the source, the CO J=6-5 emission traces the warm gas (50-300 K) 
extending to ^ 3^', and the H2 emission traces the hot gas (~ 2000 K) extending from to In 
CO J=2-l, however, the shell structure in this model is significantly less extended than that seen 
in the observations (Fig. [1]), indicating that the temperature in the shell in this model must have 
decreased less rapidly from the tip to the source than that in the observations. No HV component 
is seen in this model due to the lack of molecules at high velocity, except for the jet-like structure 
at the tip of the lobe. 

PV diagrams of CO and H2 emission cut along the outflow axis together show a V-shaped PV 
structure extending to ^ 40 km s~^ from the systemic, with the CO emission at the lower end 
and the H2 emission at the upper end (Fig. [5l model 1). This PV structure is associated with 
the LV components, arising from the shocked AGB wind in the shell, which is in the expanding 
lobe. This PV structure is expected, with the left and right parts from the back and front walls 
of the expanding lobe, respectively. In observations, the LV component in CO J=2-l also shows a 
hint of a V-shaped PV s t ructu re (Fig. [2j) and has also been modeled with an expanding lobe by 



Sanchez Contreras et al.l (120041 ) . The expansion velocity in their model, however, is ~ 22 km s ^, 
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much lower than that in our model. The upper end of the V-shaped PV structure is also expected 
to be seen in H2 even at low resolutions in the observations (Fig. [5l model 1). However, no such 
PV structure is seen in H2 observations (Fig. [J]). It is hkely because the outflow lobe in CRL 
618 indeed has a smaller expansion velocity as found in the CO J=2-l observations, so that the 
V-shaped PV structure could have been smeared out in the observations. 



4.2. Model 2: Pulsed Atomic CFW at 10,000 K 



Model 2 here corresponds to model 4 in lPaper 11 . The CFW here is the same as that in model 
1 but with a temporal variation in de nsity and velocity with an amplitude A = 0.5 (i.e., Av = 150 
km s~^) and a period P = 22 yr (see IPaper j ). In this model, a series of internal shock pairs are 
formed in the cavity moving along the outflow axis at high velocity, as the faster CFW catches 
up with the slower CFW (Figs. [St, d). They are mainly atomic. They do not affect the shell 
dynamics. They only produce small disturbances (ripples) on the shell structure in the regions 
where the internal shocks interact with the shell. Thus, the emission structures in CO and H2 are 
the same as in model 1, except that there are faint ring-like structures seen across the cavity arising 
from the ripples in the shell (Fig. [4j). The PV diagrams are also similar to those in model 1, except 
that there are faint low-velocity H2 emission near the source arising from the ripples in the shell 
(Fig. Ej). No HV component is seen from the internal shocks due to the lack of molecules. As a 
result, although this model was found to produce optical emission in the cavity with the internal 
shocks (jPaper j ). it can not produce molecular emission in the cavity. 



4.3. Model 3: Steady Molecular CFW at 1000 K 

This model is the same as model 1 but with the CFW assumed to be molecular at 1000 K, in 
order to have HV molecular gas inside the cavity. Note that the assumed value of the temperature 
is not important, because the temperature of the CFW will drop rapidly below 100 K in 200 AU 
(or 0(^2) due to radial expansion and radiative cooling (Fig. [3]f). Since the shell is not pressure- 
driven, reducing the temperature of the CFW does not change the shell dynamics (Fig. [3^). Thus, 
the molecular fraction of the shocked AGB wind in the shell is the same as in model 1, and so are 
the emission and PV structures of the LV components (Figs. [4] and [5]). On the other hand, only 
the newly shocked fast wind in the inner boundary of the shell is molecular (as indicated by the 
contours of /, Fig. [3]C) since the shock there is weak. It is hot, producing two HV shell structures 
in H2 at the far end, one at ^ 3(^5 around ~ 60 km s~^ and one at ^ 5^' around ^ 100 km s~^, 
with the velocity increasing with the distance from the source. The one at ^ 3(^5 may have a 
counterpart at ^ 4^' in the H2 observations (see Figs. [T]and[2j). Since the temperature of the CFW 
itself drops below 100 K in 200 AU, HV CO emission is seen arising from the CFW itself near the 
source. However, due to radial expansion of the CFW itself, this HV CO emission is very faint as 
compared with the LV component in the PV diagrams at the resolutions of the observations (Fig. 
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[5j), inconsistent with the observations, in which the HV CO component is comparable to the LV 
component (Fig. [2]). Besides, no HV H2 emission is seen from the CFW itself near the source. 

4.4. Model 4: Pulsed molecular CFW at 1000 K 

This model is the same as model 3 but with the CFW assumed to have a temporal variation 
in density and velocity with an amplitude A = 0.5 (i.e., Av = 150 km s~^) and a period P — 22 yr. 
As in model 2, the shell structure and dynamics are not affected much by the internal shocks (Figs. 
[3^, h). Thus, the emission and PV structures of the LV components, which arises from the shocked 
AGB wind in the shell, are also similar to those in model 2. The internal shocks, as they propagate 
beyond 1000 AU (^ 1^') from the source, become strong enough to dissociate the molecules (Fig. 
[3p). Therefore, only the internal shock closer to ~ 1" produces HV H2 emission (Figs. |4]&[5j). 
These internal shocks are so strong that the molecules in the newly shocked fast wind in the inner 
boundary of the shell at the end of the lobe are also dissociated. Thus, unlike model 3, no HV H2 
emission is seen from the newly shocked fast wind at the far end. 

4.5. Model 5: Pulsed Molecular CFW at 1000 K with A = 0.3 

This model is the same as model 4, but with a smaller variation amplitude with A=0.3 (i.e., 
/\v — 90 km s~^). In this case, the internal shocks are weaker, so that molecules can survive in 
the internal shocks except at their tips (Figs. [3, j, E, see the contours of /). Therefore, HV H2 
emission is seen arising from the internal shocks, with the intensity decreasing rapidly with the 
distance due to radial expansion (Fig. Sj). It is seen arising from the first three internal shocks, 
forming three HV emission clumps (or knots) that may correspond to the three HV H2 emission 
peaks seen in the H2 observations (Fig. [1]). The HV H2 emission seen at ^ 2" away from the source 
can be compared to that seen at similar distance in the observations. However, its kinematics, 
with the velocity increasing toward the tip and an emission peak at the highest velocity (Fig. [5j), 
is inconsistent with the observations (Fig. [2]). HV CO J=6-5 emission is seen from the wings of 
the internal shocks and may show similar PV structure to that seen in the observations. However, 
it is too weak (relative to the LV CO J=6-5 emission) to explain the observations. In addition, the 
internal shocks are too hot to have CO J=2-l emission. 

4.6. Model 6: Pulsed Molecular CFW at 1000 K with A = 0.1 

This model is the same as model 5, but with a smaller variation amplitude with A=0.1 (i.e., 
/\v — 30 km s~^). In this case, the internal shocks become too weak to dissociate any molecules in 
the internal shocks (Figs. [3j^, 1). Molecules also can survive in the newly shocked fast wind in the 
inner boundary of the shell at the far end of the lobe, producing HV H2 emission at the far end as 
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in model 3. The internal shocks are cold enough to have HV CO emission in J=6-5, but too cold 
to have HV H2 emission and too hot to have HV CO emission in J=2-l. In addition, due to radial 
expansion, the intensity of the HV CO emission in J=6-5 at ~ 2" is still too faint to be compared 
with the LV component (Fig. [5D, inconsistent with the observations (Figs. [T]and[2]). Moreover, 
there will be no atomic and ionic emission from the internal shocks, inconsistent with the optical 
observations that show atomic and ionic emission along the body of the lobe. 



4.7. Summary of our models 

In summary, limb-brightened shell structures are seen in our models in CO and H2 at low 
velocity arising from the shocked AGB wind in the shell, and can be identified as the LV components 
in the observations. However, the shell structure in CO J=2-l is significantly less extended than that 
seen in the observations. None of our models can properly reproduce the observed HV molecular 
emission near the source along the body of the lobe. In our models in which the CFWs are atomic, 
the shell has high-velocity molecular material only at the far end of the lobe. Thus, in order to 
produce HV molecular emission near the source, the CFW itself has to be molecular. In our steady 
CFW models, although there is some HV CO emission near the source from the CFW itself, it is far 
too weak (relative to the LV CO emission) to explain the observations. This is because the column 
densities of the HV molecular material decrease very rapidly due to radial expansion. In our pulsed 
CFW models, HV H2 emission can be seen along the outflow axis arising from the internal shocks 
and may correspond to that seen in the observations. However, its kinematics is inconsistent with 
the observations. HV CO 6-5 emission can be seen arising from the internal shocks and may show 
similar PV structure to that seen in the observations. However, it is too weak (relative to the LV 
CO J=6-5 emission) to explain the observations. In addition, HV CO J=2-l emission from the 
internal shocks is far too weak. 



5. Discussion 
5.1. Could the HV component be from a different lobe? 

In high-resolution optical images, CRL 618 is clearly seen with multi ple lobes on each side 



and some of these lobes are not as extended as the Wl lobe (see Fig. [T] and iTrammell fc Goodrich 



20021 ) . It is thus possible that the HV component is actually associated with a small lobe that 
happens to be aligned with the Wl lobe, and arises from the shell instead of the internal shocks. 
This small lobe could be either a highly inclined extended lobe that appears small in projection or 
a younger lobe that has a smaller physical linear extent compared to the Wl lobe. Our model 3, 
which has molecular emission from both the shocked AGB wind and shocked fast wind in the shell, 
can be used to investigate this possibility. 
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At a higher indination of 60°, the outflow lobe has a projected length of ^ At this higher 
inclination, the emission from the shocked AGB wind has a higher velocity, but still much lower 
than the observed HV component (Fig. [S]). The two HV H2 emission from the shocked fast wind 
at the far end are now projected to ^ 2'!^ at ~ 100 km s~^ and ~ 3^' at ^ 180 km s~^, respectively. 
The one at 2'!^ can be compared with the observed H2 component at the same distance. However, 
its velocity increases with the distance, inconsistent with the observations. In addition, no HV CO 
emission is seen from the shocked fast wind. Only faint HV CO emission is seen near the source 
from the CFW itself. 

At a younger age of 86 years, the outflow lobe also has a projected length of ^ 3^' at the 
inclination of 30°. The emission from the shocked AGB wind and shocked fast wind has the same 
velocity structure as at older age (Fig. [6]), because the shell is momentum-driven and thus roughly 
self similar. The HV H2 emission from the shocked fast wind at ^ 2" may show a similar PV 
structure to that seen at ^ 2'!^ in the observations. However, no HV CO emission is seen from the 
shocked fast wind. Again, only faint HV CO emission is seen near the source from the CFW itself. 

As a result, adding a small lobe into our models is still not able to reproduce the observed HV 
CO and H2 components in CRL 618 properly. The HV molecular gas in CRL 618 might have a 
layered structure, with different temperatures and thus different emissions in different layers. 



5.2. Could the CFW have a different physical structure? 



The CFW may have a different physical structure, as concluded in iPaper ll when comparing 
optical emission along the body of the lobe in our models to the observations. Here we discuss the 
two possible physical structures of the CFW that have been studied in the literature. 



5.2.1. Episodic Cylindrical jet? 



The CFW could be in a form of a cylindrical jet, with the fast wind material confined to a 
small cross section and collimated to the same direction along the outflow axis. In this case, the 
CFW has a constant density with the distance from the source. The observations in CO J=2-l 
also show that the HV compon ent is better reproduced by a cylinder with the gas flowing axially 
(ISanchez Contreras et al.ll2004l ). Simulations with episodic cylindrical jets have been performed by 
a nun iber of authors in order to reproduce the morphology and kinematics of protostellar outflows 
(e.g., ISuttner et all IiQQtI : iLee et all boOlh . Such jets can also pr oduce a series of in ternal bow 
shocks with HV CO and H2 emissions along the body of the lobe (ISuttner et al.l 119971 ). Without 
radial expansion, the emission of these bow shocks does not decrease as fast as th at in our curren t 



models. In addition, unlike those in our current models, the shocks are ballistic (Lee et al 



mull- 



The cylindrical jet models have been a pplied to PNs (ICliffe et al.l Il995l : ISteffen fc Lopea Il998l : 
Lee fc Sahai II2OO4I : iGuerrero et al.ll2008l ) and may apply to PPNs as well. 
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The jet is likely to be magnetized because the magnetic field can provide the required col- 
limation. In addition, a magnetized jet has also been foun d to better reprodu ce the jet emis- 
sion in the young PN Hen 2-90 than an unmagnetized jet fjLee fc Sahai I l2004l ). The jet could 
be launched by magneto-centrifu gal forces from a magnetiz ed accretion dis k and star system 



(IBlackman. Frank, fc WelchI I2OOII ). like the protostellar jets (jShu et alJ Il994l : iKonigl fc Pudritz 



2000). The central star could have a binary companion and th e accretion disk could form as 



the material flows to the companion (IMastrodemos fc Morrislll998l ). 



5.2.2. Bullet? 



Recently, a bullet (i.e., a massive clump) model has also been prop osed for PPN shaping 
because it can also produce a collimated outflow lob e (IDennis et al.ll2008l ). The bullet could be 
launched by an explosive MHD mechanism (see, e.g.. Matt et all 20061 ) . It is not clear, however, 
how a bullet can produce HV molecular emission near the source along the body of the lobe. One 
way that we can think of is to have a chain of bullets along the body of the lobe. Note that, 
however, a jet model with a periodic variation in velocity has been found to be better than a model 
with a chain of bullets or dense clumps, in reproducing; the kno tty jet emission in the young PN 
Hen 2-90 (e.g., comparing models 1 and 2 in lLee fc Sahai l[20o3 ). 



Dennis et al.l (120081 ). having compared their jet and bullet models, argued that the bow-shock 



heads of bullets take on a V-shaped configuration and are thus more consistent with the observations 
of CRL 618, whereas bow-shock heads of jets are more U-shaped. However, the assumed density 
and thus the mass-loss rate of the AGB wind in their models are about a factor of 1000 less than 



those assumed in our niodels , which are more appropriate for CRL 618 (IKnapp fc MorrisI Il985l : 
Sanchez Contreras et al.ll2004l ). With a high density of the AGB wind as observed, the je t model 
has been found to produce V-shaped heads (see, e.g., ISuttner et al.l 119971 : iLee et al.ll200lh similar 



to tha t seen in CRL 618, because of efficient cooling at high density. The jet model in lDennis et al 
(j2008l ) may i nstead apply to another PPN OH 231.8+4.2, which shows a wide and more U-shaped 
outflow lobe (IBujarrabal et al.ll2002l ). 



5.3. Is the CFW molecular? 



In our models, the CFW is assumed to be molecular in order to have HV molecular emission 
near the source along the body of the lobe. It could be intrinsically molecular. It could also be 
intrinsically atomic and becom e molecular rig!;ht after launched due to its high mass-loss rate, as 



proposed for protostellar wind (IGlassgold et al.lll99ll ). Alternatively, the HV molecular emission 



co uld arise from the entrained AGB material that was originally close to the source, as suggested 
bv ICox et al.l (120031). or in th e inner dense torus-like core around the source, as suggested by 
Sanchez Contreras et al.l (120041 ) . It is not clear in our simulations, however, how this entrainment 
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can happen once the sheh is formed and shields the fast wind from being interacting with the 
AGB wind around the source. It probably can happen if the fast wind is actually launched from 
around a binary companion. A s the AGB wind flows to the companion, an accretion disk can form 
(IMastrodemos fc Morrislll998l ) and launch the fast wind, and then the fast wind can entrain the 
AGB wind. 



6. Conclusions 

We have presented simulations of CFW models including molecular cooling and time dependent 
chemistry of hydrogen. We have also derived from our simulations intensity maps and position 
velocity diagrams for CO J=2-l, J=6-5, and H2 1-0 S(l) emission and compared them to recent 
observations of CRL 618. In our models, limb-brightened shell structures are seen in CO and 
H2 at low velocity arising from the shocked AGB wind in the shell, and can be identified as the 
low- velocity (LV) components in the observations. However, the shell structure in CO J=2-l is 
significantly less extended than that seen in the observations. None of our models can properly 
reproduce the observed high-velocity (HV) molecular emission near the source along the body of 
the lobe. In order to reproduce the HV molecular emission, the CFW is required to have a different 
structure. One possible CFW structure is the cylindrical jet, with the fast wind material confined 
to a small cross section and collimated to the same direction along the outflow axis. 
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Table 1. Physical parameters of the CFWs 



Model 


Temperature (K) 


Composition 


Status 




1 


10^ 


Atomic 


Steady 




2 


10^ 


Atomic 


Pulsed with A = 


0.5 


3 


10^ 


Molecular 


Steady 




4 


10^ 


Molecular 


Pulsed with A = 


0.5 


5 


10^ 


Molecular 


Pulsed with A = 


0.3 


6 


10^ 


Molecular 


Pulsed with A = 


0.1 
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Fig. 1. — CO J=2-l, CO J=6-5, and H2 1-0 S(l) emission contours superposed on the HST optical 
image of CRL 618. (Left Column) The top an d bottom panels show the con tours of CO J=2-l 
emission at high and low velocities, respectively (ISanchez Contreras et al.ll2004l ). (Middle Column) 
The top and bottom panels show t he contours of CO J=6-5 emission at high and low velocities, 
respectively (INakashima et al.ll2007l ). (Right Column) The top, middle, and bottom panels show the 
contours of H2 1 -0 S(l) emission at low velocity, middle velocities, and high velocities, respectively 
dCox et al.l[2003l ). The source is either at the position offset (0,0) or marked by a cross. The Wl 
lobe is the optical lobe extending ~ &^ to the west from the source, as indicated. 
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^LSH (km s ) 
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Fig;. 2. — Posit ion- velocity (PV) diagrams of H2 (black c ontours with g!:ray imag!:e in bo th panels, 



Cox et al.ll2003l ). CO J=2-l (red contours i n the left panel. [ Sanchez Contreras et al.ll2004l ). and CO 
J=6-5 (red contours in the right panel, Nakashima et al. 200?! ) emission cut along the outflow axis. 
The Wl lobe is mainly redshifted with a velocity ranging from —40 to 150 km s~^ LSR. 
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Fig. 3. — Simulations of our CFW models from 1 to 6, at the age of ^ 160 yrs. (a), (c), (e), (g), 
(i), and (k) show the number density of hydrogen nuclei in logarithmic scale, (b), (d), (f), (h), (j), 
and (1) show the temperature in logarithmic scale. The molecular fractions / = 0.1 (white line) 
and / = 0.49 (magenta line) are shown to outline the distribution of molecular gas. The separation 
between the AGB wind and the CFW is shown by the color tracer c = 0.5 (gray line, at which 
half is AGB wind and half is CFW). The CFW is atomic at 10,000 K in models 1 and 2, while 
is molecular at 1000 K in models 3, 4, 5, and 6. The CFW is steady in models 1 and 3, while is 
pulsed with A=0.5 in models 2 and 4, A=0.3 in model 5, and A=0.1 in model 6. (A) — (F) show 
respectively the blow-ups for the regions in the boxes A— F. 
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Fig. 3. — Continued 
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Fig. 4. — Integrated CO J=2-l (top row), CO J=6-5 (middle row), and H2 (bottom row) intensity 
maps derived from our models 1 to 6 from left column to right column. The first set and second 
set of panels show the LV and HV components at 0(^1 resolution, respectively. The third set and 
fourth set of panels show the LV and HV components at 0('5 resolution, respectively, for comparing 
to the observations. The HV components are derived by integrating the emission with velocity 
higher than 50 km s"""^ from the systemic and the LV components by integrating the emission with 
velocity within 50 km s~^ from the systemic. 
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Fig. 5. — PV diagrams of CO J=2-l (top row), CO J=6-5 (middle row), and H2 (bottom row) 
emission cut along the outflow axis, derived from our models 1 to 6. Top set is derived with a 
spatial resolution of 0['l and a velocity resolution of 4 km s"-"^. Bottom set is derived with a spatial 
resolution of 0(^5 and a velocity resolution of 5.2 km s~^ for CO and 9 km s~^ for H2 emission, 
for comparing with the observations. The contours go from 10% to 90% with a step of 20% of the 
peak value. 
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Fig. 6. — PV diagrams of CO J=2-l (top row), CO J=6-5 (middle row), and H2 (bottom row) 
emission cut along the outflow axis, derived from our model 3 at a higher inclination of 60° (left 
column) and a younger age of 86 yrs (right column). They are derived with a spatial resolution of 
0('5 and a velocity resolution of 5.2 km s~^ for CO and 9 km s~^ for H2 emission, for comparing 
with the observations. The contours go from 10% to 90% with a step of 20% of the peak value. 



